Electronic structure and optical properties of ZnS/CdS nanoheterostructures 
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The electronic and optical properties of spherical nanoheterostructures are studied within the 
semi-empirical sp 3 s* tight-binding model including the spin-orbit interaction. We use a symmetry- 
based approach previously applied to CdSe and CdTe quantum dots. The complete one-particle 
spectrum is obtained by using group-theoretical methods. The excitonic eigenstates are then de- 
duced in the configuration-interaction approach by fully taking into account the Coulomb direct and 
exchange interactions. Here we focus on ZnS/CdS, ZnS/CdS/ZnS and CdS/ZnS nanocrystals with 
particular emphasis on recently reported experimental data. The degree of carrier localization in the 
CdS well layer is analyzed as a function of its thickness. We compute the excitonic fine structure, 
i.e., the relative intensities of low-energy optical transitions. The calculated values of the absorption 
gap show a good agreement with the experimental ones. Enhanced resonant photoluminescence 
Stokes shifts are predicted. 
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The quantum confinement effects in semiconductor 
nanocrystals (NC's) or quantum dots (QD's) have been 
investigated over the last two decades. A recent develop- 
ment concerns nanoheterostructures containing concen- 
tric layers of two different semiconductors. The initial 
idea of capping a NC with a shell of higher-gap (bar- 
rier) material in order to reduce the uncontrolled sur- 
face effects has been extended to quantum-dot quantum- 
well (QDQW) structures with an inner layer of lower-gap 
(well) material which allow a greater control over the op- 
tical properties. [l| There are several examples such as the 
inclusion of a layer of HgS in a CdS NC [lj, |2j,|3j, where 
the lattice parameters of both compounds are of similar 
size (a = 5.851 A for the /3-HgS and a = 5.818A for 
the CdS). Recently, ZnS/CdS based QDQW's have been 
synthesized 0, 01 with a larger lattice mismatch (more 
than 7 % , the ZnS lattice constant being a — 5.409). 
Some other combinations have also been studied, such as 
CdSe/ZnS or CdSe/CdS.@ 

On the theoretical side, effective mass approximation 
(EM A) models were proposed 0,0 without a full many- 
body treatment of the Coulomb direct and exchange in- 
teractions. An atomistic theory was needed, however, 
to adequately describe a QDQW where the well can be 
as thin as one monolayer. Very recently we proposed a 
tight-binding (TB) model [jj, which is a generalization 
of our previous work. This theory successfully accounted 
for the optical properties of CdSe [HJ and CdTe [H03 
NC's. Independently, Bryant and collaborators have also 
developed a TB approach. 0, 0] In contrast with the 
Bryant group, we include the spin-orbit interaction, take 
into account the symmetries and obtain the full single- 
particle spectrum before deducing the exciton states. 
The Coulomb and exchange interactions were already in- 
corporated in our study of CdS/HgS/CdS QDQW'sQ, 
prior to the more recent work by Xie et aL[lJ] To our 
knowledge, the present paper reports the first TB cal- 
culation of the low-energy excitonic states in CdS/ZnS 
nanoheterostructures. After a brief description of the 
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theoretical model which can be completed with our pre- 
vious publications [flt llflll2l |. we present our main results 
on ZnS/CdS, ZnS/CdS/ZnS and CdS/ZnS nanocrystals, 
including a comparison with the available experimental 
data. 

Our calculations are based on the semi-empirical TB 
model for bulk semiconductors introduced by Vogl et al. 
[l5| and generalized by Kobayashi et al. to account for 
the spin-orbit interaction. [ 
orbitals are used: s, p x , p y , p z , s 
ical s* orbital was introduced in order to simulate the 
effects of d orbitals on the conduction bands. The in- 
teratomic hopping matrix elements are restricted to the 
nearest neighbors. The model is then described by a 
total of 15 independent parameters. The CdS TB pa- 
rameters used here are essentially those proposed earlier 
by Lippens and Lannoo 0] which we have slightly mod- 
ified in order to incorporate the spin-orbit coupli ng. The 
ZnS parameters are those calculated by Bertho.[lj| As 
for the band offset, following Ref. the valence band 
maximum of ZnS is set to eV and that of CdS shifted 
to +0.4 eV. The bulk gaps of ZnS and CdS are 3.7 eV 
and 2.5 eV respectively. 

We outline the method for deriving the symmetrized 
TB Hamiltonian (see Ref. E for details). The NC's, of 
roughly spherical shape, are constructed starting from 
a cation at the origin by successively adding nearest- 
neighbor atoms through tetrahedral bonding. We pas- 
sivate the NC surface by placing a hydrogen s orbital 
at each empty nearest-neighbor site on the surface dan- 
gling bond so that the final surface states are several eV 
far from the gap edges. Hereafter we use a simplified 
notation: We write r co /r w /r c i to indicate that a NC is 
built with a core radius r co , with r w and r e i designating 
the widths of the well and clad layers, respectively. The 
lengths are all given in angstroms. 

The single-particle TB Hamiltonian is reduced to a 
block diagonal form by writing it in a symmetrized ba- 
sis corresponding to the double-valued representations 
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Tk (k — 6,7,8) of the Td point group. flfl| This method 
not only reduces the size of the Hamiltonian to diago- 
nalize, but also provides the exact symmetry classifica- 
tion of the wavefunctions. One can thus deduce selection 
rules and related insights over the relative intensities of 
optical transitions between the valence and conduction 
states. 0| We then diagonalize the Hamiltonian and ob- 
tain the full one-particle spectrum in a given nanostruc- 
ture: all the energy levels and eigenstates. We have also 
calculated the spatial projections of the density of states 
(DOS) onto the core, well and clad regions, when re- 
quired, to study the influence of the position and thick- 
ness of different layers on the whole energy spectrum. 

The optical properties concern exciton-like elemen- 
tary excitations which can be described in terms of the 
Coulomb direct and exchange interactions between the 
electron and hole. The total Hamiltonian can be written 
as, [13 

Hvc^v' c' — £-v)&w' &cc' Jvc.v'c' K vcv ' c ', (1) 

where e c , — e v are the electron and hole energies, J and K 
are the Coulomb and exchange interactions respectively 
(see Ref. E3 for more details). Note that we leave the e— h 
exchange interaction (K) unscreened up to the nearest- 
neighbor site (primitive cell of the zincblende crystal), 
in analogy with the unscreened short-range part in the 
theory of bulk exciton. 

The unscreened on-site Coulomb and exchange inte- 
grals for anion (cation) are assumed to be U cou i — 20(6.5) 
eV and U exc h — 1(0.5) eV respectively. These values fol- 
low roughly those obtained for CdSe and CdTe. [13, E3 
The values for the S atom in CdS have been calculated by 
means of a simple scaling. We assume the same parame- 
ters for ZnS. The nearest-neighbor exchange integrals are 
assumed one tenth of the on-site ones. The permittivity 
is taken as e(oo): e(oo) = 5.7 for ZnS and e(oo) = 5.2 
for CdS in single binary QD's. We take an average value 
in ZnS/CdS, e(oo) = 5.4. Strictly speaking, the dielec- 
tric constant in a nanocrystal is size-dependent and sub- 
stantially reduced with respect to the bulk low-frequency 
value e(0). See, for example, Ref. |2lj. where a nseudono- 
tential calculation for CdSe nanocrystals is presented: 
The dielectric constant increases with increasing size and 
asymptotically approaches the bulk high-frequency value 
e(oo). To our knowledge, no such calculations are avail- 
able for the present nanoheterostructures nor even for 
their components. We, therefore, assume an average of 
the bulk high-frequency values. 

The large lattice mismatch in the ZnS/CdS system is 
expected to lead to elastically strained structures with 
the relatively thin shell/well layers conforming to the 
substrate lattice constant. Q Unfortunately, there is no 
simple way of including this effect in our TB model. Pre- 
liminary calculations suggest that a simple modification 
of the interatomic matrix elements alone, following the 
Harrison scaling rule for the bond length, is inadequate: 
It yields a systematic increase of the optical gap in the 
ZnS/CdS and ZnS/CdS/ZnS structures considered here 



by 10-30 %. The effect is, of course, qualitatively inverse 
in CdS/ZnS. A more refined modelization scheme for 
structural distortion involving the diagonal matrix ele- 
ments (see, for example, Ref.[22]) seems necessary. Work 
is in progress in this direction and will be reported else- 
where. The results presented below have been obtained 
by neglecting the strain effects. 

In Fig. we show the total density of states for some 
ZnS, CdS and ZnS/CdS/ZnS NC's. It can be seen how 
the projection of DOS onto the CdS well quickly grows 
when one and two monolayers are added. A monolayer 
yields a CdS well projection, almost equivalent in size 
and energy distribution, to the core ZnS projection. The 
addition of a second layer leads to a CdS well projec- 
tion which is much larger than the corresponding ZnS 
projection. Interestingly, the size of the gap is governed 
essentially by the ZnS properties. We also look at the 
radial distribution of charge corresponding to the eigen- 
states. We can see the numerical values of the probabil- 
ity of presence of electron and hole in the LUMO and 
HOMO in TableQ]for several cases. The idea behind the 
comparison is to investigate the role of the successively 
added CdS layers. Both the electron and hole show an en- 
hancement of their presence in the well region; the effect 
is more pronounced for the hole. When two CdS layers 
are added the hole is essentially trapped in the well. In 
Table ITT1 we show the two highest valence (HOMO) and 
two lowest conduction (LUMO) levels. The ZnS single- 
particle spectrum presents an intrinsic degeneracy which 
is the consequence of the absence of spin-orbit coupling: 
When a level is Tj — Tg degenerate its orbital symme- 
try is T§, because the spin- 1/2 representation is Tq and 
T$ x Tq = Tr + T$. Similarly, the degeneracy — T$ 
indicates that the underlying orbital symmetry is T4. 

The excitonic spectrum is calculated by diagonalizing 
the exciton Hamiltonian in Eq. (1) in the configuration- 
interaction method with as many valence and conduction 
states as necessary to reach numerical convergence. The 
relative intensities of optical transitions are next com- 
puted by following the procedure described in Ref. (TrT |. 

In order to discuss the available experimental data, let 
us first note that in QDQW structures like CdS/ZnS: 
in Ref. Q the dispersion for dots with a mean diam- 
eter of 58 A is a — 9 A. On the other hand, in small 
nanocrystals, theoretical results show big jumps in the 
energies when one adds just one or two shells to a given 
NC: See Table ITT1 Our procedure to fix the adequate 
size for a given set of experimental data 0] is as follows: 
We first analyze several NC's of sizes differing by one or 
two atomic shells around the experimental size and retain 
three such structures yielding a good agreement with the 
experimental values of size and optical absorption gap. 
The procedure is illustrated in Fig. [21 We show the fine 
structure of absorption for three different ZnS NC's of 
around 20 A diameter. We can see that the experimen- 
tal onset from JJ], indicated by a vertical dashed line, is 
in good agreement with our results if we remember that 
the onset must have a finite width. Note that the succes- 
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sive sizes in Fig. correspond to the addition of single 
atomic shell. In Fig. we show the fine structure of 
a series of QDQW's based on a ZnS core of radius 7.2 
A. The experimental values are given in the figure cap- 
tion. Finally, it is worth mentioning that the predicted 
photoluminescence Stokes shift is enhanced by the pres- 
ence of the barriers. In particular, the Stokes shift in the 
R = 14.8 A CdS NC is 12 meV, whereas the 7.2/5.9/2.6 
QDQW, even slightly bigger in size, shows a shift of 53 
meV. 

In Fig. 0we show a comparison between ZnS and CdS 
cores with ZnS/CdS and CdS/ZnS core-shell structures. 
The results are in accord with the experimental data of 
Ref. Q: The absorbance changes little between CdS 
and CdS/ZnS in the low-energy region, but the change 
is drastic between ZnS and ZnS/CdS. The former result 
has to do with a surface that is already fully passivated 
so that a thin barrier shell enhances the localization only 
slightly. On the other hand, in ZnS/CdS, the outer shell 
is a well and the spatial distributions are drastically mod- 
ified, leading to a different absorption spectrum. 



To conclude, we have studied the electronic structure 
and optical properties of spherical nanoheterostructures 
based on ZnS and CdS within a symmetry-based tight- 
binding approach. The full single-particle spectrum, the 
low-energy excitonic states as well as the relative intensi- 
ties of optical transitions are calculated. A careful analy- 
sis of the spatial distribution of the DOS and the HOMO 
and LUMO wavefunctions reveals that only two mono- 
layers of CdS suffice to trap both the electron and the 
hole in the well. Our principal results concern the exci- 
tonic fine structure yielding the absorption gap and the 
resonant photoluminescence Stokes shift. Their size and 
composition dependences have been analyzed on a mono- 
layer scale. We obtain a satisfactory agreement with 
the available experimental values of the absorption on- 
set without using any adjustable parameter. The pre- 
dicted Stokes shifts indicate a strong enhancement in the 
nanoheterostructures as compared to simple nanocrys- 
tals. 

This work has been supported by the Spanish Minis- 
terio de Ciencia y Tecnologla project MAT2002-00699. 
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Tco 1 T-w 1 T c l 


Core 


Well 


Clad 


7.2/3.2/0.0 


0.592(0.323) 


0.383(0.671 


0.0(0,0) 


7.2/5.9/0.0 


0.322(0.118) 


0.659(0.879) 


0.0(0.0) 


7.2/2.7/2.9 


0.419(0.171) 


0.328(0.590) 


0.241(0.237) 


7.2/5.9/2.6 


0.211(0.046) 


0.623(0.663) 


0.158(0.289) 



Table I: The probability of presence of electron(hole) in the 
LUMO(HOMO). The core radius and the well and clad widths 
are given in the first column for each case (in angstroms). 



r co /r w /r c i (A) 


HI (eV) 


H2 (eV) 


LI (eV) 


L2 (eV) 


10.1/0.0/0.0 


-0.322(8, 7) 


-0.378(8,6) 


4.251(6) 


4.598(7, 8) 


7.2/3.2/0.0 


-0.156(8) 


-0.167(6) 


3.995(6) 


4.291(8) 


7.2/5.9/0.0 


0.056(8) 


0.036(6) 


3.662(6) 


3.922(8) 


7.2/2.7/2.9 


-0.046(8) 


-0.058(6) 


3.844(6) 


4.112(8) 


7.2/5.9/2.6 


0.188(8) 


0.166(6) 


3.574(6) 


3.809(8) 


0.0/14.8/0.0 


-0.247(8) 


-0.256(8) 


2.988(6) 


3.335(8) 



Table II: ZnS/CdS/ZnS: Energies of the two highest (lowest) 
occupied (unoccupied) molecular orbitals, denoted by HI, H2 
(LI, L2). The symmetry index n designating T n is shown in 
the parenthesis. The size and composition are given in the 
first column. 
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Figure 1: Density of States for three NC's. Top: DOS in a ZnS QD of radius R = 10.1 A. Middle(bottom): A ZnS/CdS 
nanocrystal with a core of R = 7.2 A and one (two) CdS monolayer (s). In the middle and bottom figures the shaded region 
corresponds to the DOS projected onto ZnS core region and the solid line to the DOS projected onto the CdS well region. 



Figure 2: Excitonic fine structure in ZnS QD's. The experimental absorption onset is indicated by the vertical dashed line. We 
show three similar sized QD's of around 10 A in radius. 



Figure 3: Excitonic fine structure in QDQW's. The core radii and shell thickness are all given in A. The respective experimental 
absorption onsets are: (a) 3.85 eV, (b) 3.55 eV, (c) 3.76 eV, (d) 3.45eV. 



Figure 4: Excitonic fine structure in core and core-shell nanostructures. The core radii and shell thickness are all given in A. 



